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We show that the boost variable, the conjugate to the coordinate rapidity, which is associated 
with the center-of-mass motion, encodes the information about the finite size of colliding nuclei in 
a Lorentz-invariant way. The quasi-elastic forward color-changing scattering between the quantum 
boost states rapidly grows with the total energy of the collision and leads to an active breakdown 
of the color coherence at the earliest moments of the collision. The possible physical implications of 
this result are discussed. 
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> ' I. INTRODUCTION 

o ' 

It is commonly accepted that on the scale of the strong interaction, which is responsible for nuclear integrity 
and compactness, the large nuclei have a macroscopically finite size and a well-defined boundary.^ This size can be 
physically measured in the rest frame of a nucleus, and it undergoes the Lorentz contraction in the moving frame 
' without any physical limitations (as is required by special relativity). In this paper, we suggest to take this fact as 

in ' 

ON 

ultrarelativistic energies. Of these consequences, the most important is the change of the symmetry: The incoming 

in 
o 



a guideline, and explore the consequences of the finite size of the nuclei for the quantum process of their collision at 



nuclei are prepared in a homogeneous space having a given energy and momentum. The fixed space- time point of 
the first interaction corrupts the initial symmetry, and enforces a different choice of the conserved quantum numbers 
""■-^ ' for the later stages. Of the ten symmetries of the Poincare group, only rotation around the collision z-axis, Lorentz 
O ,' transformation along it, and the translations in the transverse x and y-directions survive. Therefore, it is profitable 
to choose, in advance, the set of normal modes which have the symmetry of the localized initial interaction and 
carry quantum numbers adequate to this symmetry. These quantum numbers are the transverse components, pt, of 
momentum and the boost, v — p°z — p z t, of the particle (which is associated with the center-of-mass motion and 
' replaces the component p z of its momentum). 

These geometric considerations can be reinforced by the quantum mechanical arguments. Indeed, from the per- 



spective of an external observer, the first thing that happens during the collision is a precise measurement, by means 
of the strong interactions, of the collision coordinate within a very short time interval. Therefore, statistically, by 
the uncertainty principle, the secondaries with any conceivable momentum p z can be detected after collision. This 
is a well-known scheme of the Heisenberg microscope. The higher resolution we want to achieve, the larger must be 
the energy resources of the microscope. In the textbook example of the electron probed by the photon, the electron 
receives energy from the hard photon. In nuclear collisions, both the kinetic energy of the nuclei, and the energy of 
the compression of the Lorentz-contracted nuclei are used for the purpose of a precise measurement of the coordinate. 
An internal observer that penetrates the future of the collision with the nuclei will see a violently expanding matter 
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x By the finite size, we mean the size which is measured by means of the strong interaction of two nuclei. If the primary 
interaction were electromagnetic (as is in the ep-or eA-processes) , then the whole concept of a finite size would become 
doubtful. 
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around him. The two viewpoints perfectly complement each other. The short scales of primary interaction provide a 
sufficient motivation to use the wedge dynamics which describes the fields inside the future domain of the "wedge" 
t 2 = t 2 — z 2 > 0, and employs the "proper time" r as a Hamiltonian time of the evolution and the coordinate rapid- 
ity 77 as a longitudinal coordinate [^JQ|. The infamous rapidity plateau persistently observed in high-energy nuclear 
collisions strongly supports this picture. 

The approach advocated in this paper explicitly incorporates the macroscopic finite size of the interacting objects 
into the quantum theory of the earliest stage of the collision. We assume that there is no measurable gluon fields 
outside the large stable nuclei. Consequently, the time moment and the z coordinate, along the collision axis, of the 
first interaction are defined with the accuracy of at least ~ 0.01/to, which is both the size of a Lorentz- contracted 
individual nucleon, and the characteristic scale of color correlation in the z direction before the collision. The full 
size of the Lorentz-contracted gold nucleus at the energy ~ lOOGeV per nucleon is ~ 0.1/m. We show, that despite 
an almost infinite Lorentz contraction and the quantum nature of the interaction process, the information about 
the finite size of the incoming nuclei does not fade away. It remains clearly identifiable in terms of the properly 
chosen Lorentz-invariant variable, the boost, which is associated with the center-of-mass motion. Thus, it is possible 
to describe the collision of the two nuclei staying on the same physical ground in any reference frame, either in the 
reference frame of one of the nuclei, or in the laboratory frame where both nuclei move almost at the speed of light. 

The fact that nuclei have finite size is intimately connected with the gauge nature of the strong interactions. 
Therefore, when addressing the problem of interaction of the two compact nuclei, we must refer to the properties of 
the vector gauge fields .[^ The colored sources of these fields must be located inside the nuclei and they can be physically 
resolved only after the two nuclei overlap. This is the only assumption we make regarding the internal structure of 
a nucleus. By all means, location of a material object inside a nucleus implies that its center-of-mass should move 
with this nucleus without crossing its boundary. Therefore, before the collision, it is natural to characterize such an 
object by its center- of-mass, i.e., by its boost v. The valence quarks are the first candidates to be considered in this 
manner. In this sense, we follow the idea of McLerran-Venugopalan model J3| in the form given by Kovchegov and 
Mueller 0]. However, we do not try to populate the nuclei with the wee partons. We think that they are gradually 
created in the course of collision |^|,^| . 

The framework of wedge dynamics also offers a unique opportunity to avoid various technical problems encountered 
when the moving at the speed of light nuclei, V z = ±c, are taken as the first approximation |^,^|. This state cannot 
be reached as a continuous limit of V — > c and a significant effort has been made to smooth out the singular behavior 
of quantum fields at V = c ^-^|. The wedge form of Hamiltonian dynamics is free of this difficulty. Furthermore, 
the gauge A T = of the wedge dynamics can be fixed completely. Hence, the transverse and longitudinal fields are 
well separated and the gluon propagators of wedge dynamics have no spurious poles that can stimulate a singular 



2 Addressing the issue of interaction of finite-size nuclei, one should keep in mind the source of the major difference between 
QED and QCD phenomena. The local gauge symmetry of QED can be extended to a global gauge symmetry which generates 
the conserved gauge- invariant global quantum number, the electric charge, which can be sensed at a distance. The proper field 
of an electric charge is the main obstacle for the definition of its size. On the other hand, the radiation field of QED appears 
as a result of the changes in the extended proper fields of accelerated charges, and one can physically create such an object as 
a front of electromagnetic wave. In QCD, the local gauge invariance of the color group cannot be extended to is global version 
that would correspond to a gauge- invariant conserved charge. Hence, we can readily define the size of the colorless nucleus, 
but we cannot create a front of color radiation in the gauge-invariant vacuum. These two properties of QCD both work for us. 
They allow one to use the Lorenz contraction to localize the primary domain of the collision and thus, to impose the classical 
boundary conditions on the color fields at later times. The existence of the collective propagating quark and gluon modes at 
the later times is the conjecture that has to be verified by the study of heavy ion collisions. 
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behavior of scattering amplidudes g. In this framework, one can use the same dynamics and the same gauge for the 
description of both incoming nuclei and the products of their reaction ||[| , thus avoiding all glitches of the "on-line" 
changing the gauge and redefinition the states ||]. 

Below, we concentrate on a specific interaction in the expanding system that emerges in ultrarelativistic nuclear 
collisions. It is mediated by the longitudinal part of the gluon field. It seems to be the leading one at the earliest 
moments of the collision of the two nuclei, and to result in the intensive color exchanges even in quasi- elastic 
subprocesses. Eventually, these exchanges must cause an active breakdown of the fragile color coherence of the 
colliding nuclei and stimulate intense color radiation. The rate o\ of these color exchanges between the quantum 
boost states appears to be large at the earliest moments of the collision, and it grows as log 2 E with the total energy 
E. This major result of this paper is given by Eq.( 3.27 ). The log 2 E dependence of the rate on the total energy of the 
collision resembles the obtained in the early sixties estimate on the maximal rate at which the total cross section may 
grow with the energy. It is known as the Froissart bound, and a steady growth of the total cross section is indeed 
observed in the proton-proton collisions. 

Originally, the Froissart bound was derived in the scope of the axiomatic field theory, a powerful approach based 
on the most general requirements, like Lorentz invariance, causality, unitarity, completeness of the basis of physical 
states and the cluster decomposition principle (see Refs. |hJ for the details). Since the perturbation theory (usually in 
a given order) can lead to an anomalously large total cross section (and thus to apparently violate unitarity) it is said 
that the perturbative total cros-section requires unitarization. Recently, this problem received a vigorous attention 
in connection with the evolution equations for large nuclei at low x p ]ll]J| . A physical protection from an excessive 
growth of cross section due to collinear problems was offered in Ref. H. From this standpoint, one can infer that the 



result (3.27) of this paper indeed complies with the unitarity. Though this issue has to be studied in more details, we 
suggest a plausible simple physical arguments below. 

The axioms of unitarity and completeness clearly are not truly ndependent. Discussion of any issue related to 
unitarity requires that the spaces of the initial and final states are completely specified. Physically, this means that 
the measurement is not accomplished until its products are analyzed. What the particular states are, depends on 
the detectors that resolve these states. In nuclear collisions, one cannot rely on the conventional "external" distant 
detectors. The role of the detectors for the earliest subprocesses (which only very tentatively can be viewed as the 
independent acts of scattering) is played by the subsequent interactions. The next-to-best thing one can do is to try 
to answer the following question. Let the fields excited at the beginning of the collision are expanded over a system of 
states characterized by some quantum numbers. Let two such states interact. What is the rate of these interactions? 
The answer will be related to the two main problems of ultrarelativistic heavy ion collisions. First, the known rate of 
the primary interaction will help to estimate the entropy production. At this point, the explicit knowledge of the final 
states is imperative, because the entropy is the number of the excited degrees of freedom. Second, it will be directly 
connected to the total cross section. Indeed, if the fields change their colors during the time ~ 1/E with sufficient 
probability, then the nucleons will lose their coherence and fall apart. A new composition of hadrons will be created 
with the probability one, and it does not really matter how the interacting states are chosen. This argument has been 



3 The division of the gauge field into the longitudinal and transverse parts can be done only with respect to the property of 
propagation: transverse fields are emitted and then propagate being limited in space-time by the light- cone boundaries, while 
the longitudinal fields are simultaneous (in terms of the Hamiltonian time) with their sources. In QCD, this scheme can be 
practically implemented only in the framework of perturbation theory, which is assumed throughout this paper. 
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tested long ago: the total cross section of the e + e~-annihilation into hadrons coincides with the cross section of the 
process e + e~ — * qq. One of the recently studied examples is the interaction of the eikonalized quarks or gluons fl2| . 
In this paper, for the same purpose, we consider the "natural" states of the wedge dynamics, deliberately leaving 
the key question of what interacts at the very beginning of the collision open. We find that, because the states of 
wedge dynamics carry internal currents in the coordinate rapidity direction, there exists a specific contact interaction 
of these currents which grows when r — > 0, and leads to the amplitude of interaction, proportional to log E. [The 
contact term in the gluon propagator has been singled out in the course of the complete fixing of the gauge A T = 0, 
and its main effect is confined to the nearest vicinity of the light wedge, r = 0, where the boundary conditions that 
fix the gauge are imposed.] If the QCD indeed falls under a jurisdiction of the axiomatic field theory (which by no 
means is self- evident), then our perturbative result, which exactly reaches the Froissart bound, may point to the 
major physical mechanism that triggers the scenario of ultrarelativistic heavy ion collisions. 

The paper is organized as follows. In section |l| we introduce the variables of wedge dynamics and clarify the 



physical meaning of the boost in classical and quantum contexts. In section III we use the boost states to estimate 
the amplitude of forward scattering with color transfer at the earliest moments of the collision, paying attention 
to the contact interaction in the expanding system. In Appendix A, we demonstrate that the contact term has 
no counterparts, and that the standard Coulomb-type terms are still there in the propagator. They are somewhat 
modified, just in a way which one could expect on purely physical grounds. In Appendix B we show, that the 
contribution of the other terms into the forward scattering amplitude is subleading. 

II. CLASSICAL AND QUANTUM PARTICLES IN WEDGE DYNAMICS 

In this section, we address the basic connection between the classical and quantum aspects of the interaction of 
compact relativistic objects, in order to prepare the stage for a more involved analysis of the interaction picture. 
First, we discuss the role of the classical Lorentz boost as a natural variable which, by its origin, is closely related to 
the finite size. Second, we review the meaning of the boost as a quantum number, and establish its connection with 
the classical boost. Finally, we show that the genuinely classical distribution of the boosts in stable nuclei before the 
collision plays a role as the initial data for the primary quantum interactions between nuclei. 

A. Introducing the variables. 

The wedge form of relativistic dynamics works inside the future domain of the hyperplane t — z = (light wedge) 
were the two finite-size ultrarelativistic objects touch each other for the first time. The natural coordinates inside 
this domain are parameterized by the proper time r and the rapidity coordinate rj, 

t = t cosh??, z = t sinhfy . (2-1) 

In terms of these variables, the action for a classical particle is 



S = J Ldr=-m J ds = -to J y/l - v 2 dr = -to / dr\J 1 - r 2 ?) 2 - r , (2.2) 
where v 2 = T 2 rj 2 + is the spatial velocity squared, and the dot means derivative over the (Hamiltonian) time r.[j] 



'Following a tradition, we use the Greek indices for the four-dimensional vectors and tensors in the curvilinear coordinates {rj 
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The canonical momenta of this particle, 

dL mT 2 f] _ dL mr 

9?/ y/l - V 2 dr \/l - v 2 

are conserved by virtue of the equations of motion. The Hamiltonian is of a standard relativistic form, 

fn 

H = vfi+p-r-L= -== , (2.4) 
V 1 — v 

which, after excluding the velocities, can be rewritten in terms of the canonical momenta, 



H = ^m 2 +P 2 + (2.5) 

The useful relations of geometric origin, which will be often referred to later on, are 

1 / sinh?7 cosh 77 3 \ m t . , . 
p 1 = -— p v = - 1 — — p — p° I = — — smh(?j - 9), 

H = cosh 77 p° — SU1I177 p 3 = m t cosh(r] — 6), (2-6) 
where m 2 = m 2 + p 2 , p° = m t cosh#, and p 3 = m t sinh# are the Cartesian momenta. Therefore, the boost, 

v = p v = ttth smh(r) -6)= x 3 p° - x°p 3 = p°(z ~ V z t), (2.7) 



is related to the center-of-mass coordinate. According to Eq.(2.5), the quantity vjr plays a role as a local longitudinal 
momentum. 

The Hamilton- Jacobi equation for the classical action of a particle reads as 



dS 1 fdS\ 2 fdS^ 2 



dr y t 2 \dri J \dr 
It allows for the separation of variables and has a solution 



+ m 2 = 0. (2.8) 



S = vti + p '■ f — I \j m 2 H — -dr = vri + p ■ f — \ m 2 T 2 + v 2 + vArsh . (2-9) 

In a quantum context, this action serves as the phase of a semi-classical wave function, tp ~ e zS , with the quantum 
numbers v and p, either when v ^> rrrit or when Tint > v. An isolated solution with the not separated variables is 

S =p- f- m t T cosh(?7 - 9) . (2.10) 

It corresponds to a plane wave, and its parameter, the (momentum) rapidity 8, is not a canonical momentum. 



is an exception, it always stands for the rapidity direction), and the Latin indices from a to d for the vectors in flat Minkowsky 
coordinates. We use Latin indices from rtom for the transverse x- and y-components (r, W = 1, 2), and the arrows over the 
letters to denote the two-dimensional vectors, e.g., k = (k x ,k y ), \k\ — kt. The Latin indices from i to n (i,...,n = 1,2,3) will 
be used for the three-dimensional internal coordinates u l = (x,y,rf) on the hyper-surface r = const. 
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B. Classical trajectories. The physical meaning of the boost v. 



In order to understand the physical meaning of the boost variable v, the canonical conjugate to the rapidity 77, 
one has to figure out how it enters the classical equations of motion. According to the Jacobi theorem, the action 
S(x n ,a n ), known as a function of coordinates x n and arbitrary constants a n , allows one to find an additional set 
of the conserved quantities, dS/da n = b„ . While the constants a n — dS/dx n usually are the canonical momenta 



corresponding to the cyclic coordinates and conserved due to the equations of motion, as in Eq. (2.9), the constants 



b n appear to be the initial coordinates. Applying the Jacobi theorem to the action (2.9), and choosing the constants 
in such a way, that at r = we have x = xq, and that at r — > 00 we have rj — 9, we obtain the equation of the particle 
trajectory, 



x(t) - X 



Px 



7](t) 



T 2 m; + v 1 — \v\ 



-Arsh- 



(2-11) 



Despite their unusual appearance, these two equations parameterize a straight line, as it should be for the free motion 



of a point-like particle. Let us rewrite the second of equations (2.11) in two ways, 



zp 



tp z 



and 



mtT sinner) — 9] 



m t T cosher) — 9] = tp — zp z = \ T 2 mf + v 



m t z* , 



m t t* , 



(2-12) 



(2-13) 



where the arrows point to the special choice of the reference frame with 9 — O.f] Then the first of the equations ( 2.11 ) 
becomes 



x(t) 




— -> — (U - \z*\) 

mf m t J nit 



(2.14) 



obviously satisfying the required boundary condition at r = 0. Now, it is easy to understand that the quantity 
vjm% is the r- independent coordinate 2* of the particle in the co-moving frame. By the definition, this quantity 



is Lorentz- invariant: the boost v is the same in all Lorentz frames. The Cartesian form ( 2.14 ) of the trajectory is 
obviously continued to all quadrants of the tz plane. This classical definition of the boost is fairly operational but, 
as the reader may notice, it requires that the base world line (plane), from which the coordinate z*, is measured is 
explicitly chosen. For the two colliding nuclei, it is natural that the base lines (corresponding to the rapidities ±Y) 
go through the point t = z = 0, where the nuclei touch each other by their surfaces. If the nuclei have radius R 
and are built from the fragments of the (transverse) mass m t , then the boosts for the right-moving nucleus will be 
in the range —2mtR < v < 0, and in the range < v < 2mtR for the left-moving one. There is no contradiction 
with quantum mechanics at this point, since the nuclei are macroscopic stable objects which can be kept under 



5 We consider the physical design of the nucleus as almost static, and neglect the possible velocity V* of the nuclear constituent 
in the nuclear rest frame. In any case, it cannot be large without undermining the alleged stability of the nucleus. The origin 
of the transverse mass may be different. It includes both the Lagrangian mass and the "adjoint mass" due to the transverse 
momentum. Inside a stable nucleus, the momenta most probably characterize the standing waves which are not likely to be 
too short, if the nucleus is in the ground state. 
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non-destructive control (in their co-moving reference frames) before the collision. Asymptotically, they have the well 
defined rapidities 9 — ±Y, which can be also measured classically, without any contradiction with the anticipated 
uncertainty relation, AvAr/ > 1. Indeed, the boosts v ss m t z* are measured inside the nuclei, while the measurements 
of the velocities of the nuclei is performed by external devices. Therefore, the boost variable is indeed perfectly suited 
for the description of the finite size objects. If the relative boosts of all constituents do not change in the course of 
the interaction, then the object remains unaltered in its (possibly new) rest frame. 




FIG. 1. Geometry of a nucleus-nucleus collision in the center-of-mass reference frame (left) and in the rest frame of one of 
the nuclei (right) The dark gray lines correspond to a semi-classical boost state in the right-moving nucleus before the collision. 



As the matter of fact, the boosts provide an invariant measure of the distribution of the constituents of the compact 
objects. The picture of rectilinear trajectories holds outside the light wedge also. Therefore, the classically prepared 
distribution of the boosts is resolved as the distribution of the further interacting quantum states with the given 
boosts, when two such objects collide. Though Eq. (2.12) expresses the boost v via the invariant mt and distance z*, 



in a quantum picture, the boost v is an independent conserved additive quantum number. 

For isolated point-like (and thus, structure-less) objects, the practical measurement of the boost requires that the 
rapidity r\(r) is measured at two time moments along the same trajectory. Then, solving the system of two equations 
( [2.12 ), one finds v and 9, the boost and the asymptotic rapidity of the particle. It is unrealistic to perform such 
measurements with sufficient accuracy in the asymptotic domain of the macroscopically large r. Unlike the case of 
the macroscopic finite size object, this kind of measurement does meet quantum-mechanical obstacles. 



C. The boost v in quantum context. 



The quantum-mechanical measurement of the boost v is very similar to the measurement of a usual momentum 
and relies on the definition of the operator of the boost, 

. d 



drj 



(2.15) 



as the operator of translations in the 77-direction. Then it becomes evident, that a simultaneous measurement of 
coordinate rj and momentum v is limited by the uncertainty relation, 



Av An > 1 



(2.16) 



In the field-theory formulation, the boost operator is given by the generator of the Lorentz rotations in the tz plane. 
In the internal geometry of wedge dynamics, the boost operator is given by the TV) component of the energy- momentum 
tensor. The boost of the quantum field at the proper time r, 



T T, \x)Tdr]d 2 r= / d£ M M^ 03 (x) 



(2.17) 
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(where M^ v \ — x v T^\ — x\T^ v + is the usual angular momentum tensor) is the integral of motion corresponding 
to the translation symmetry (Lorentz rotation) in r\ direction. The quantum states with the given boost v are the 



eigenstates of the operator ( 2.15 ), and their eigenfunctions depend on n as e lvri . The full wave function of a scalar 
particle with the boost v and the transverse momentum p is the solution of Klein-Gordon equation with the separated 
variables r, r\ and ft, 

*$( x ) = ^^H^l(m t r)e^ + ^ . (2.18) 
It is normalized on the hypersurfaces r = const, 

a 



i>;>,?{x) i>8,p(x) Td V d 2 r = 5(9 - 6')S(p-0) . (2.19) 

zonst 

This equation normalizes the measurements performed by an array of the detectors moving with all possible velocities. 
At any particular time of the Lorentz observer, this array even does not cover the whole space. 

At large v ^> 1, and v > m t r, which is relevant to the earliest stage, the asymptotic of this solution is semi-classical, 

„27r/4 ~ivr]+ipr 

42 (*) « l-T Kr2 + , 2]1/4 e-V^+^ArshW^) K e iS (2 2Q) 

clearly indicating that at the small time r the quantum particle with the finite boost v continues to follow its classical 
trajectory, since its classical action is large. Indeed, the surface of the light wedge, everywhere except for its vertex, 
corresponds to n — > 00. 



The wave functions (2.18) are connected, by means of Fourier transform, with the plane-wave solutions, 



The saddle point of the Fourier transform (|2.2l| ) (or its inverse) is located at the value of v (or 9) defined by the 



relation, v = Tm t sinh(# — n), corresponding to the classical definition (2.7) of the boost. One can easily see that 



these wave functions also are semi-classical with the action (2.10), and have a usual momentum (or the rapidity 9) as 
a quantum number. These states become localized in rapidity rj at later times, Tm t S> 1, and namely these states are 
most likely to be detected by the expanding collective system. 

The key element of the suggested approach is that the Lorentz-invariant boost states, which are independently 
prepared in the two approaching nuclei, begin to interact as the quantum states only when the nuclei overlap. At 
this moment, the positions of the nuclei constituents (classical boosts, which describe the elementary constituents of 
the nuclei even outside the light wedge) are translated into the quantum numbers, which define the periodicity of the 
wave functions in the coordinate rapidity direction. It is evident, that at the earliest times, the distortion of the initial 
geometric picture should be only minimal. Therefore, it will be a sufficient approximation to study the transitions 
into other boost states, and we stay within this approximation until the end of this paper. The rate at which these 
early distortions develop appears to be quite large. 

The dynamics of boost states preserves the invariant information about the finite size of the nuclei both in the 
laboratory frame when each of the two nuclei is contracted up to a negligible small size, and in the rest frame of one of 
the nuclei (target) when the second one (projectile) passes through it as a seemingly infinitely sharp shock front. One 
cannot assign to the moving in the x + direction front a finite width, neither in the z-, nor in the ^"-direction, without 
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a conflict with the special relativity. On the other hand, in the framework of the wedge dynamics which operates with 
the boost states, it is safe to consider the limit of the infinite momentum frame at the end of the calculations. H 



III. SCATTERING IN WEDGE DYNAMICS 

The nuclei meet each other at the two-dimensional plane t = z = 0, where the first interaction take place. This 
interaction resolves the nuclei constituents (e.g., the "partons", or "color dipoles") with the boost v as 0, and excite 
the quantum states with the boost v ~ 0. The wave functions of these states do not depend on the rapidity coordinate 
?7, and evenly fill in the interior of the light wedge. At the same time, the two precursors, which are most likely to 
be the fronts of the propagating gluon field, begin their way in the light-like directions, t ± z = 0, thus creating 
the physical boundaries of the light wedge, t 2 = t 2 — z 2 = 0. Passing through the nuclei, the precursors resolve 
the elements with the finite boosts, which are negative for the right-moving nucleus and positive for the left-moving 
one, and initiate a transient process of interaction between the nuclei. These interactions excite the quantum states 
with positive and negative boosts, which depend on r\ as e wn . In this way, the classical boosts, v c i = mtZ*, are 
transformed into the quantum numbers of the wave functions which have the period 2it/is in the 77 direction, and 
occupy the entire future domain of the point t = z = 0. Before the collision, the nuclei as a whole, are the coherent 
states of QCD, and their (color) coherence cannot be destroyed immediately. At r — > +0, the resolved boost states 
have the same phases they had in the nuclei: the decomposition of the nuclei in terms of the boost states is still a 
coherent superposition. ^Furthermore, since the classical action of the states with the finite boosts is large, even the 
resolved partons continue to move along their rectilinear classical trajectories. The character of the further evolution 
crucially depends on the subsequent interactions. Below, we study the quasi-elastic forward scattering of the colored 
quarks prepared and detected in the given boost states. This scattering is mediated by the gluon field and results in 
the color exchange which alone is capable of destroying the coherence of the nuclear wave function. 

The propagators of the gauge fields in wedge dynamics were studied in Q,^]. In Appendix A, we review their 
properties with the emphasis on the needs of the present study. The leading contribution comes from the spatially 
local "contact term" of the longitudinal part of the propagator. In order to give a flavor of its origin, we have to 
emphasize, that we study the phenomenon where the finite charge density is formed as a result of the interaction, 
and the proper fields of the gradually created and yet delocalized charges physically overlap with their sources. Thus, 
aiming at the dynamic picture, we have to give the priority to the currents, expressing the charge density p(t) via 
the divergence of the current, dtp = — V • j, which eventually generates the contact term in the propagator. The 
effect of the evolving charge density p(t) becomes fully included into the Hamiltonian Hint = j • A, which is the only 
form compatible with the completely fixed gauge A T = 0. This evolution of the color charge density is the result of 
the interference between various partial waves, and it is not connected with the motion of the physically resolved 
point-like color charges. Without an interaction, these partial waves would coherently sum and form the stable nuclei. 
Of those interactions that take place when the nuclei intersect, the most important are the ones which lead to the 



This, however, leaves open the question of what is detected in high-energy collision. The answer crucially depends on what 
kind of the quantum mechanical ensemble is involved in a particular measurement. 

7 The boundary condition A v {t — 0) = imposed on the gauge fields in the wedge dynamics, together with the gauge condition 
A T — 0, makes it impossible that the fields of precursors immediately modify the phases (rotate the color charges) along the 
light-like planes x + = and x~ = 0. This property, which allows one to "switch on" the interaction between the nuclei without 
an artificial color-changing "shock wave," is in contrast with the case of the null-plane dynamics with the gauges A = 0. 
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largest transition amplitudes. 

An apparent complexity of the formulae in the wedge dynamics is caused by the curvature of the hyper-surfaces 
of the constant r. The hyper-surface r = +0 is the one where the initial data are naturally set, and in has an 
infinite curvature. An explicit dependence of the internal metric on t makes the vector differential operators more 
cumbersome and leads to an interplay between the longitudinal and transverse fields. 



A. Choosing the observable 



Wedge dynamics deals only with the fields that emerged from the localized collision of two macroscopic objects. 
This collision is considered as a precise measurement of the partons coordinates at the finite time moment r — ► +0. 
Therefore, it is impossible to pose a formal scattering problem with the asymptotic initial states. Instead, we take an 
approach based on the calculation of the Heisenberg observable 13 j| , 



N(l', 2') = (1, 2|n(l')(n(2') - S vv )\l, 2) = (0 1 a 2 ai S^, 4 a v a 2 /Sal 4 10) , (3.1) 

which is the inclusively measured number of pairs of the final state field excitations with quantum numbers 1' = k[) 
and 2' = (i^k^) (fc includes the transverse momentum and boost, i - color). This observable is evolved from the 
initial state of the two interacting field excitations with quantum numbers 1 = and 2 = (i2>fe)- This quantity 

is closely related to the total cross section. Indeed, assume that the measurementis an impulse process that freezes 
decomposition of the colliding nuclei in terms of the eigenfunctions of the corresponding operator. This decomposition 
can become incoherent only due to real interaction, which will either excite the new states, or just break the phase 
balance between the initial ones. All this will contribute to the probability that the initial state is altered, i.e., to 
the imaginary part of the forward scattering amplitude. The color exchanges take place at the earliest possible time 
Tmin ~ 1/ y/s, and create a new color composition which must eventually (with the probability one) evolve into a 
new composition of hadrons. We emphasize, that a particular choice of the basis of the interacting at r > boost 
states is important only as long as we are interested in the rate at which the color coherence is broken. The color 
transfer between the boost states seems to be extremely intensive at the beginning of the collision. The geometry of 
the collective modes which will be the actual final states can be quite different p4| . 



Expression (3.1) is bi- linear with respect to the evolution operator S and thus it cannot be treated according to the 
Feynman rules. For its evaluation on should use the so-called Schwinger-Keldysh technique Jlq ] in the form adjusted 
for the calculation of inclusive amplitudes uM . The evolution operator for the problem of evolution of the observable 



(3.1) is of a usual form, 

S = Texp{i J H int {x)d A x) (3.2) 

with the Hamiltonian 

H int (x)=f(x)A fl (x)=f(x)[AW(x) + f dzDlT 3] {x,z)f{z)\ , (3.3) 



where the second term in brackets is the longitudinal field A^° n9 \x). The propagator D%,° n9 \x, z) implicitly contains 
6(x° — z ). For the sake of defmiteness, consider the fermion color current, 

j»(x) a =g* i (x)tl jl »* j (x) , (3.4) 

and commute the final-state Fock operators with S and using the commutators, 
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di(k)S — Scii(k) = I dzip^^z) -= 



St a J(fc)-o+(fc)5t = | rfo^^^)^) , (3.5) 

In this equation, (z) is the one-particle wave function from the decomposition of the field operator, 

= Y,Mk)4 +) (x) + bt(k)4-\x)) . (3.6) 
k 

These commutations result in (disconnected pieces are omitted) 

N(l',2') = [ dxidx 2 dyidy2^\y2)4> ( J\yi)(0\a2ai 



5%, {yi)5^Ay2 



a\a\\Q)4f , {x 2 )^\x l ) . (3.7) 



Here, the functional derivatives over ^ act from the left, and the derivatives over ^ act from the right. Next, we 
compute the functional derivatives retaining the terms up to the order g 2 . This yields, 

N(l',2') = g 4 / dx 1 dx 2 dy 1 dy 2 4>\f > {y2)4k/\yi) (0| a l2 (k 2 )a h (fci) 



xT^{[-%Mr%Axxh v Af^{x 2 )AW^ 
xT{^ j2 (y 2 )^^ yi )YA^(y 2 )A^ 

+i^ j MDlT a]h ' a \yuy2)i v *M] } 4(*i)4(fca) |o) ^\x^{ Xl ) . (3.8) 

The calculations are accomplished as follows. The fermion operators are contracted with the remaining Fock 
operators of the initial state, producing the final-state wave functions, and making the final adjustment of the color 
indices. This can be done in two ways, which differ by a full interchange of the quantum numbers of the one-particle 
initial states. The vacuum average of the products of the transverse gluon field operators gives the transverse part 
of the T-ordered propagator I?! 00 ! (y 2 , y\) and of the Tt- ordered propagator Dl 11 ! (x 2 , x\). | The two terms, with 
D^ on9 \y 2 ,yi) and D^ on9 \yi 7 y 2 ) cover two complementary domains, y 2 > y? an d y 2 < Vii respectively. Together, 
they form the longitudinal part of the T-ordered propagator D^°\y 2 , y\). Finally, the transition probability can be 
cast in the form, 



N(l',2')=g 4 



d Xl dx 2 l4V (z 2 faW* (sifr'Vg (si) 



x^W^iK^*?^ " the same(fci,«i <-> k 2 ,i 2 )} 



(3.9) 



8 In this paper, we use the Keldysh-Schwinger formalism |15j in its modified form developed earlier with the view of application 
to the inclusive and transient processes. We employ the notation used in Refs. . The indices of the field correlators with 

the Keldysh contour ordering of the field operators (like D[ab]) as well as the labels of their linear combinations (like Dr re ti) 
are placed in square brackets. 
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B. Scattering of scalar quarks with the given boosts v. 



The observable number of couples, N(V, 2'), can be rewritten by introducing the full set of the intermediate states 



into the Eq. (3.1 



N(l',2') = ^2(0\a 2 aiS^al,al,\X)(X\ava 2 'Sa[al\0) = ^ \{X\a v a 2 ,Sa\a 2 \Q)\ 2 . (3.10) 
x x 

In the lowest order of the perturbation theory, there is no additional emissions, and only the vacuum state \X) = |0) 

contributes, 



N(l',2') = |M 1)2 ^, 2 f 



(3.11) 



In the lowest order, the inclusive transition probability ( |3.10|) is just the squared modulus of the matrix element 
depicted on Fig. 0. 





(a) (b) 
FIG. 2. The forward (a) and backward (b) amplitudes of the qq scattering. 



1. Scattering amplitude 



Consider the matrix element of the scattering amplitude 

= g 2 J dx^^Xi^^^D^xx,^) . (3.12) 
where we exchange the spinor quarks for the scalar ones, and accordingly replace 

using the states of scalar quarks with the quantum numbers k — (fc, v), transverse momentum and boost. In this case, 
the wave functions are of the form, 

,-W2 



4>) = W 



V'fc, (z) 



2 5 / 2 7T 



H£'{mtT)e 



—iwq—ikr 



Using the propagator in the mixed representation, 



Df:\x u x 2 ) 



(2ir) 3 ™ ' 2 >^^J C 



and integrating over the spatial coordinates, we obtain 

9 2 



Mi.a-n'.a' = -%r +v2-v'i- + fc 2 ~ K - k' 2 ) 

poo poo 

Jo Jo 1 3 

Xff H (ri) 5 "" n (7a)(* 1 +k' 1 )i{k 2 + k' 2 ) m D [ Z ] {r l ,T 2 ^q) , 



(3.13) 



(3.14) 



(3.15) 
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where £ = v\ — i/[ = v 2 — i/' 2 , q = k\ — k[ = k 2 — k' 2 , and we introduced 3-vectors, g u (r)pi — {—p, —v/t 2 ) , as well as 
a short-hand notation, m 2 = m 2 + k 2 , m[ 2 = m 2 + k 2 . 

Computing the transition amplitude ( |3.12 ), we will be interested in the states with the large boosts \v\ ^ 1, 
v > m t r, \ui — V2\ 3> 1 ■ In this case, the asymptotic of the Hankel functions reads as 



^fl^H = [ne-™! 2 H\l\m t T)} 



TTU/2 TTWl 



'2m 



[ m 2 T 2 +Jy 2]l/4 



-iyj m^T 2 +^ 2 +it , Arsh(^/mtr) 



(3.16) 



We have mentioned already, that in this limit, we have ip^\x) °c exp(iS' c /), where S c i is the classical action, found 



in Sec. |l| In the limit of \ 3> rriiT, we have 



i/Arsh = v In 

mT 



\v\ ln(2|^|) — //ln(mr) , 



and the product of the four Hankel functions in the integrand of Eq. ( 3.1 5| ) becomes 

4( TOl )-^ (m 2 )-^ K)»j (m' 2 )< i(|i/l|+|l/2| _ |v ; 



xe 



TT 2 \^ 2 ^ 2 \^ 2 

-i[|i/ 1 |ln(2|i/i|)+|i/ 2 |ln(2| ! / 2 |)-|i/J|l l i(2|^|)-|i/2|ln(2|^|)] ( ^1 

T2 



(3.17) 



The last factor here is the most significant for future analysis. The rest is just the phase factor. 

In what follows, we compute the leading term corresponding to the contact part of the gluon propagator (see 
Eq. (|A.26|) in Appendix A), 



l -S(r))S(f) . 



(3.18) 



It is local in rj and r, and the modulus accounts for both terms with D^ long ^ in Eq.(3.8). In this approximation, the 



matrix element (3.15) becomes 

^1,2->1', 2' 



2 — 5(u t +v 2 -v[- v' 2 )5(hi + k 2 - k[ - + + gia 7 



(3.19) 



2(2tt) 3 ' " z ' L 'a/-v-x , - -i .-2/ 4 |j /li/2l/ /mi/2 

where a is an inessential real phase. In the approximation given by the equation ( |3 . 1 7 ) it absorbs all the dependence 
on the transverse momenta. Now, it remains to compute the integral 

-iC 



I = 



poo pOO 

/ ndn T 2 dT 2 g^(T 1 )g^(T 2 )[D^\r ll T 2 )] i 
Jo Jo 



dn / drM 



n 



2riT 2 V r 2 



(3.20) 



where the cutoffs are introduced in order to isolate the possible singular behavior. Computation is straightforward, 



21 



rdr 



TCI 

T 2 



dx[x iC + x~ i<: 



1 



The cutoffs r„ 



C 2 + 4 
tq and r, 



7^ sii.Klnir/r,,)] 



ji 2 



C 



2(l-^(l + cos[Cln(T/r )]) } , 



(3.21) 



T in this formula are the external physical input. Making a choice for tq and T 



it is useful to keep in mind that the interaction ( 3.18 ) is due to a non-stationary part of the longitudinal (Coulomb) 
field of the charges resolved at r = 0. Similar cutoffs are needed in a stationary part. A proper choice leads to the 
Coulomb logarithm in the collision term in the QED plasma, and we follow this example. 
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The only field A n that contributes the contact term (3.18) vanishes at r = 0, and its effect on the resolved at 



t = charges cannot be instantaneous. Therefore, the lower limit tq is related to the earliest time when the boost 
states belonging to the incoming nuclei are resolved by means of the strong interaction. Practically, this is the time 
which it takes two nuclei to completely overlap. Therefore, this minimal time is defined by the velocities of the 



incoming nuclei in the laboratory frame, tq ~ At this time, the stationary phase of partial waves (2.21) 



u>^\j(x), corresponding to the particles with the given rapidities 9, are stretched over the widest rapidity interval 
Ar/ ~ 2ln(y/s/m±). This estimate coincides with the well known kinematically allowed width 2Y of the rapidity 
plateau, 2Y ss log(s/m% har ).(See Ref [jjj for the further details.) 

The upper limit T has to be set because at some time T m ax the picture of the independent collisions breaks up. The 
"final" state fields are not emitted into the free space any more (which affects even the QCD evolution equations 0). 
Therefore, T corresponds to the time when subsequent interactions begin to erase the memory about the origin of 
the boost states from the compact nuclei. By this time, the system must develop collective interactions which result 
in the effective masses of the plasmon-like modes in a dense medium. It is clear that these masses can emerge only 
gradually HQ]. An attempt to evaluate this gradual process in the scope of wedge dynamics has been undertaken 
in Refs. |^,@. This calculation relies on the following physical mechanism: The low-p t mode of the radiation field 
acquires a finite effective mass as a result of its forward scattering on the strongly localized (and formed earlier) 
particles with q t ^> pt- Regardless of what the exact value of this "screening mass" is, it seems reasonable to take 
T ~ 1/ (j>d, which is consistent with the semi-classical approximation, Tfio -C v. 



The two limits of the Eq.(3.21) are of special interest. Let y/s — > oo, while ( is kept finite. Then 

T 2 



c 2 



2*5(0 -I}, (3.22) 



the amplitude is strongly confined near the forward region and the corresponding cross section diverges. 

Next, let us consider the physical limit of the forward scattering, £ — > 0, while keeping y/s finite. In this case, we 
have 

/~T 2 [ln(iyi) - 1] , (3.23) 
the inclusive amplitude is proportional to the maximal width of the rapidity plateau, Y oc hi(y/s), which is the only 



geometric factor that can accompany the contact interaction fl3.1S| ) . Its square naturally sets the upper bound for the 
scattering probability. 

The second term in the forward scattering amplitude ( |3.9| ) , which corresponds to the complete exchange of the two 
initial states (backward scattering), is obviously small. Indeed, this case corresponds to v\ ~ v' 2l and v 2 ~ v[. In this 
case, 

|C| = \vx-v[\ « \v x -v 2 \ » 1 , 



and the function ( |3.2l| ) is small. 

2. Scattering probability 

Since we consider the processes which develop in the course of a single collision, the notion of the cross section is 
not well defined. In order to deal with the quantity which is as close as possible to the standard cross section, let us 
introduce the "normalization volume" fl = 7ri? 2 K, the product of the transverse area and the length of the rapidity 
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interval over which the nuclei become expanded by the first measurement of the collision coordinates. The wave 
functions of all states begin to occupy this volume when the two nuclei have completely overlapped, i.e., by the time 
Tmin ~ l/y/s. The wave functions ip^ , given by Eq. ( 3.13| ), in the matrix element (3.9) thus acquire an additional 
factor (27r) 3 / 2 Sl -1 / 2 . The quantity p = SI -1 will play the same role as the flux factor j = 1/ ST = v re i/V in the case 
of the standard 2 — > n scattering (see, e.g. Ref. Multiplying the squared modulus of the matrix element (3. It) 

by the densities of the final states, £ld 2 k' dv' / ^it) 3 / 2 , and replacing one of the delta functions by f2/(27r) 3 ,we arrive 
at the differential inclusive probability, 



dw 



5(vi +l>2 — v'x — v' 2 )8{ki + k 2 — k[ — k 2 ) 



{v x + v' x Y(v 2 + v' 2 ) 



I \2 



Pd^dv'^k'zdv'z 



(3.24) 



2tt 16\i>iu 2 u[p 2 \ 

Dividing dw by the density p = Q -1 , we obtain the closest analog of the cross section which can be introduced in 
order to characterize a single event, 



du i = 8{u\ + 1/2 



fcj_ k 2/ 



\V2 



2n 



(3.25) 



16\vii^v[i> 2 \ 

Since I 2 has the dimension [length] 4 , the quantity a\ also has the dimension of area. The upper limit T max — T in 
Eqs. (|3^20|)-(^23|) can be estimated from the condition t/id ~ 1 < and is related to the formation of the (final) 
states as they are detected by the subsequent interactions at the later period of the evolution. In the limit of a nearly 
forward scattering, and integrating d 2 k' 2 dv 2 with the aid of the delta- functions, we arrive at 

2 



d<Ji 
d 2 k' t d( 



2tt 9 



{2 Vl + Q 2 {2v 2 ~ty 



1 



161^2(^1 + 0(^2-01 Ml, VC 2 +4 
sin[Cln( v / i//i D )] 1 + cos[( \n(^fs / p D )\ 



C 

In the limit of the forward scattering it becomes 

da i 



d 2 q t d( 



at 2 
8tt 9 



1 , 2 Vs 

— In 



(3.26) 



(3.27) 



where qt ~ k' t is considered as the transverse momentum transfer. Our basic approximation implies that this transfer 
is small, qt < pd- The color trace 



1 1 
3 ' 3 



accounts for the processes with and without color transfer. 



IV. SUMMARY 



The main result of this note is given by the Eq. ( 3.27 ). The logarithmic character of the answer (the color-changing 
amplitude cx a s ln(r m i„/r maa; ) « a s \n(y/s/ pn) ) is due to the dimensionless-ness of the rapidity and the boost 
variables, rather than to the Coulomb nature of the interaction. This answer indicates, that we may expect a massive 
breakdown of the color balance in the colliding nuclei at the earliest time r ~ 1/ ^/s. The rate at which the intensity 
of this breakdown grows with the energy is proportional to In 2 s. 
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The key assumption that led to this result is the existence of a sharp boundary of the colliding nuclei. If this 
assumption is not correct, then there is no any reason to consider the problem of the nuclear collision in the framework 
of wedge dynamics, and the whole picture of the collision will look differently This would also undermine alternative 
approaches to the problem, like the McLerran-Venugopalan model SB. An immediate logical consequence of the 
finite size is the absence, inside the stable nuclei, of the finite color charge density, which could significantly fluctuate 
and produce the long-range fields. Only under this assumption could we safely discard the static component of the 
gauge field which would correspond to the finite charge density at t = and to consider the creation of color charges 
in the course of the nuclear collision as a transient process. The currents in rapidity direction, which we relied upon 
in our calculations, appear as a result of the phase shifts in the system of delocalized fields (and thus propagating 
with the phase velocity), rather than due to the motion of the resolved point-like charges. 

The wedge dynamics was conceived as a tool which is adequate to the earliest stage of the collision, where the 
initial color coherence becomes broken. It is not applicable to the ep— DIS, where the electron probes the long range 
electromagnetic fluctuations in the proton ||,||. In its turn, the evolution equations which describe QCD fluctuations 
that accompany ep-DIS do not seem to be relevant to the collisions of the finite-size nuclei. The primary breakdown 
of color coherence in a nuclear collision (in terms of the states of the wedge dynamics, it is indeed the earliest process) 
must result in color radiation which can exist only for a short period (in proper time), only before the fields begin to 
build up the collective modes of the expanding continuous media. Only these collective effects can bring in the scale 
(/i£>) to the entire process and serve as a feed-back that limit the intensity of the primary emissions PJl^]. Later on, 
the dynamics of the process must become local on this scale. 

The transient process of the plasma formation will come to its saturation at the moment when the growing with 
time (and density) effective masses of the collective modes begin to screen all emission, from the evolving sources, at 
the scales below the one given by the dynamically generated effective masses ^,0- Being unable to radiate, these 
sources must pass through and form the receding nuclear remnants. Thus, it is likely that the total energy of the 
collision is responsible only for the time-scale of the initial interaction and the full width of the rapidity plateau, while 
the parameters of the final state in the central rapidity region are universal and independent of the initial energy 
(above a certain threshold). Eventually, the total energy of the colliding nuclei is shared by the newly born matter 
and these receding remnants. It is not clear yet if the quark-gluon matter will have time to sufficiently thermalizc 
and be described by a single parameter, the temperature. However, it seems unavoidable that the entropy created 
at the earliest moments must result in the pressure, which is the first thing we shall try to theoretically estimate. A 
success at this point will much simplify the whole scenario by allowing to incorporate the hydrodynamic picture from 
a sufficiently early proper time. 

Our preliminary estimates show that the boost states of wedge dynamics do not effectively scatter with large 
transverse momentum transfer. Further analysis is necessary to verify this estimate, which (being correct) could 
explain the absence of high-p t jets observed in the first available RHIC data. The jets are not strongly quenched, they 
well can be absent at all! Does this mean that perturbative QCD is totally unrelated to the ultrarelativistic heavy 
ion collisions? We do not think so. It just has to be used in a different way than in ep-DIS or pp-collisions. The 
major source of this difference has been first outlined in Ref. ||: in nuclear collisions, the final states that saturate the 
unitary cut in the ladders that correspond to QCD evolution equations cannot be saturated by quark and gluon states 
in free space. In this paper, we point to the fact that the initial states can be different from the free massless wee 
partons with the given light-cone momenta. They well can be the boost states of valence quarks which are explicitly 
confined inside the finite-size nuclei before the collision. 
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V. APPENDIX A. THE GLUON PROPAGATOR. 



In order to study the interaction of the two charged states with the given boosts v, one needs to have an explicit 
form of the gauge field propagator. Particularly, since the boosts are additive and obey the conservation law, one 
need to know what the quanta of radiation that carry the boost quantum numbers are. It is also necessary to know 
the form of the proper (longitudinal) fields produced by the charged particles. In this section, we present a detail 
analysis of the gluon propagator in wedge dynamics, which has been derived in Ref. The main purpose is to 
carefully trace the origin of the new contact term. At first glance, it may look abnormal since it neither shows up 
in the field of a moving static charge, nor has it any properties associated with the propagation. We want to show 
that all Coulomb-type terms still exist in the propagator. They are somewhat modified, in a way which one could 
expect on purely physical grounds. Namely, the Coulomb fields vanish outside the future domain of the point where 
the charge was created. Our analysis indicates that the other parts of the propagator cannot hide anything similar 



to the exclusive contact part which is solely responsible for the final result, Eq. ( 3.27 ), of this paper 



A. The field of a static source 



The field of a static source in wedge dynamics is found || when one solves the linearized (Maxwell) equations 
of motion without the external current, imposing the gauge condition A T — 0. An additional boundary condition, 
which allows one to fix the gauge completely, is A^t = 0) = 0. In fact, this condition brings nothing new, since 
the hypersurface r = is light-like, and the r- and ^-directions are degenerate there. In this way, one finds three 
modes, of which two, V^ E \x) and V^ M \x) are the transverse fields. The modes V^^ and V^™> are normalized 

vk vk 

according to a usual definition of the scalar product in the functional space of the solutions of the Maxwell equations, 

(V,W)= [ d v f d 2 fTg lk V*id T W k , (A.l) 

and satisfy the Gauss law without the charge. The third mode, y( stat ) has zero norm, and its definition is accomplished 
with the aid of Gauss law with the static source. [In the absence of any currents, the source can be only static] The 
electric and magnetic fields of this mode are 

dvd 2 k e ^v+ikf 



Ef tat \r,r,ri) 



(2tt) 3 



ik 2 



vk 2 TS-i,iv{k t T) 



p(k,u) 



(A.2) 



B\ stat] (T,r,r,) 



dvd 2 k e ivr ' +ikf 
(2^) 3 k\ 



hy 

~k x 




us-i tiv (ktT)p(k,u) , 



(A.3) 



where s m ,iu{x) is the Lommel function, a solution of the inhomogeneous Bessel equation with x m 1 as the external 
source, 

f" + -f'+(l + ^)f = x m - 1 . 
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There exists an extremely important relation between the two Lommel functions,! 

s\,i v {k t T) + v 2 S-i <iu (k t r) = 1 . 



(A.5) 



First of all, it is necessary in order to verify that the electric field of a static charge distribution, ( A.2), indeed satisfies 
Gauss law,R 



-d v E v + rd r E r = rj T = p 



(A.6) 



Second, it is precisely the unit on the right side of Eq. (A.5) which will give rise to the contact term in the full 
propagator. 



The Fourier component of the vector potential of the static field is 



k r Q-\,iv{ktT) 

(27r) 3 ^ t 2 [ vQi,i v {hr) 



(A.7) 



where we introduced the functions 



(*e) 



,{x) dx 



In spite of an unusual appearance, this is nothing else but Coulomb's law in the framework of wedge dynamics. In 
order to see this explicitly, let us consider the system of point-like charges located at the points fi in the transverse 
plane and moving with rapidities 9^ 



rj T = ^qi S(n - 9i)5(r- n). 



(A. 



For a single charge, the explicit form of the electric field components is 



Ef ati (T,f,r,) 



q_ 

47T 



? cosh(r; — 9) 
2 sinh(?7 — 6) 



Q(r - r t ) 

R 3 



f/r 2 
tanh(?; — 9) 



S(r-r t ) , 



(A.9) 



where R = [r 2 + r 2 sinh 2 (77 — 9)] 1 ^ 2 , is the distance between the points (0, 9) and (r, rf) in the internal geometry of 
the surface r = const. On can obtain the first term in this formula taking the usual (gauge-independent) expression 
for the electric field of the moving charge, transforming it to the new coordinates, and multiplying it by the 9(r — rt), 
which eliminates the field outside the light cone of the point where the charge had emerged. The second term (with 
the light cone delta-function) corresponds to the wave front that accompanies the process of the charge creation at 
r = 0.[] 

Since the electric field is Ei — d T Ai, the vector potential is recovered by means of integration, 



9 It is useful to keep in mind the integral representation 

si iu(ktr) = 1 / cos (ktT sin cj>) cosh i/(f>d(/> , (A. 4) 

shin -kv 

which indicates that the functions s\^ v (x) and v 2 ! S-i l i v (x) are regular at v — 0. 

10 In terms of the physical components, £ m = ^/—gg mn E n = %/ —gg mn d T A n , the Coulomb law reads exactly as in Cartesian 
coordinates, d m £ rn = p. 

n This is not a true radiation. The real Coulomb mode A\ stat ^ is orthogonal to the complex propagating modes vf TE ' > and 



(TM) 
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Ai(t)= / E^dr' / ^(r')rfr'. (A.10) 



o 

Now, when r t is taken as the actual lower limit, the result of the integration explicitly coincides with the Fourier 
transform of Eq. ( |A.7| ). The Fourier transform of the Lommel functions appears to be discontinuous in an exactly 
relativistic way (the details of its calculation are in the next section) . 

One may ask how the Coulomb mode could be found from Maxwell's equations of motion which do not include 
Gauss' law. The answer is simple and natural: the Coulomb field outside the static charge distribution must satisfy 
the equations of motion for a free field. 

There are two surprises connected with the static solutions of the wedge dynamics. First, the source is static if 
it expands in such a way that its physical component J T — tj t (t, 77, r) does not depend on t. Indeed, the charge 
conservation has its physical form, d^J^ 1 = 0, only in terms of the physical components = s/—gg^ v j v of the 



electric current. The second surprise is the light-cone boundary of the static field in Eq. (A. 9). 

Finally, let us consider the conservation of the charge of a fundamental field in full QCD. Now, the equation of 
charge conservation reads as 

d^ + gf abc AlJ? = 0. (A.ll) 

Let only the j T a component of the current differs from zero. Then for the charge Q a = J Tdnd 2 rj T , we have 

d T Qa + gf abc A b T Q c = . (A.12) 

Since the gauge condition is A T = 0, we conclude that Q a = const. In the framework of wedge dynamics, the notion 
of static charge is well defined even if the individual charges move with respect to each other (in a specific way). 
Similar result can be obtained in the system with the Hamiltonian time t — x° with the gauge condition A = 0. If 
all (color) charges are at rest, their proper static field does not "rotate" their color. However, this will not be the case 
if we chose a different gauge condition, e.g. A 3 = or divA = 0, which would require that A ^ 0. 

Finally, it is easy to understand, that since the proper gluon field of the static fundamental color charge does not 
affect the charge itself, this gluon field cannot be a carrier of the color charge. This fully agrees with the fact, that the 
norm of the Coulomb mode equals zero, because its field is real (contrary to the complex fields of the transverse modes 
that represent gluons). An additional reason to pay special attention to the static field configuration is that the field 
corresponding to the charge density p(r = 0) is an isolated exceptional static field. It was necessary to describe it in 
details in order to have a reference point for a more involved analysis of the fields created by the charged currents. 



B. The full longitudinal field 



The gluon propagator, which we review and analyze in some details below, was found as a (retarded) response 
function between the potential and the current for the linearized (Maxwell) equations of motion. The potential is 
represented as a sum of three terms, 

A = A^ + A^ + A^ inst ^ = A^ + A^ ™ 9 ^ . 

The second and the third terms constitute the longitudinal (in a sense of the Gauss' law) field. The goal of this 
somewhat technical analysis is to demonstrate that the longitudinal part of this propagator indeed includes a new 
contact term. At the same time we want to show, that the standard Coulomb fields are still present in the propagator, 
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almost unchanged and are modified only by the relativistic causal boundaries which one would expect to appear for 
the fields of the emerging charges. 

The transverse part of the retarded propagator is trivial. It is built from the partial solutions of the homogeneous 
wave equations, 



A? r \ Xl ) = / d 4 x 2 e( Tl - r 2 )A' i t ; i ) (x 1 ,x 2 )j m (x 2 ) 



where 



dv d 2 k V \V^ ) (x)V^ ) *(y)-V^ ) *(x)V^ ) (y)} 



(A.13) 



(A.14) 



X=TE,TM 



which can be easily recognized as the Riemann function of the original homogeneous hyperbolic system. The Riemann 
function solves the boundary value problem for the evolution of the free radiation field. It is obtained immediately as 
a bilinear expansion over the full set of solutions of the homogeneous system. 
The name of the instantaneous part is motivated by its explicit form, 



Af nst \k,u; Tl )- Pik '^ n) 



(2n) 3 ik 2 



k r Q-i,iu(ktTi) 
vQi,iu(hn) 



(A.15) 



the 'potential A 1 ." 1 ^ is simultaneous with the charge density p = rj T . Formally, it can be obtained by adding the time 



dependence to the charge density in the expression for the static potential (A. 7). However, this form is inconvenient as 
long as we have to use A a j^ = Aij 1 as the basic form of the interaction Hamiltonian. Therefore, we have to eliminate 
the charge density p completely, and replace it by the spatial components j n of the current. The replacement follows 
an evident prescription, 



^ dp 

p(ri,v,k) - p(0,v,k) = I dr 2 — 

o ot 2 



T 2 dr 2 [k s j s (r 2 , v, k) + v.f(T 2 , v, k)] 



(A.16) 



The effect of the initial charge density po = p(0, v, k) would correspond to the clearly visible static pattern in the 
longitudinal part of the field. In the framework of perturbative QCD, this pattern is not active, since it cannot 
transmit the color charge. Furthermore, as we have argued previously, in nuclear collisions, the initial density of the 
color charges at r = is zero. This leads to 

rTl r 2 dr 2 



Af nst \ Tl ,v,k) = - 



In this form, the three remaining (in the gauge A 
via the spatial components of the current. 

The dynamical longitudinal field A^ is of the form 

" TI r 2 dr, 



k r Q-i^ v (k t Ti) 
vQ\,iv{ktT\) 





k s 


l 


V 



j m {T 2 ,V,k) 



(A.17) 



0) spatial components Ai of the vector potential are expressed 



Af\r l ,v 1 k) = 



kf 
V 



k s Q-is v {k t T 2 ) 

vQl,iv{k t T 2 ) 



3 m {T2,v,k) 



(A.18) 



It also does not allow for the bilinear expansion with two temporal arguments. Its electric and magnetic field is 
simultaneous with the current j m also. In what follows immediately, we intend to single out the contact part of the 
propagator, which shows up only in the D nv component and connects A v with j v . 

In order to set the stage, it is instructive to start with the electric and magnetic fields of these two modes, E m =A m , 
£ m = V~ gg mn An and B m — —(2y / —g)~ 1 e mln Fi n . Since the potential A\ L ^ is the three-dimensional gradient, we 
immediately see that BY^ — 0. [Note, that A^ is the gradient of a time- dependent function, and thus is not a pure 
gauge.] Starting from the expression for A^ lnst \ and using the relation [^J, 
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-l,iv(faT) ~ Ql,iv(k t T) = -^^-Si <iu (k t r) = T-^-8-n v (ktT) 

V 1 OT OT 



we obtain by a straightforward calculation that 

rT T 2 dT 2 V 

o (2^» 



B? nst] (r,v, k) 



ky 




k x 




k x 




ky 


S-l,i v {k t T 







V 


m 



(A.19) 



(A.20) 



i.e., the longitudinal part of the magnetic field has only the azimuthal component (the magnetic field circulates around 
the current flowing in the ^-direction) , which is natural for the distribution of charges which experience expansion in 
z-direction. Note, that the magnetic field exists even when p is r-independent. 
In the same way, we compute the electric fields 



e\ L] (t,v, k) 



{2-Kfkl 







V 


I 



k B Q-i,iv(hr) 

vQ\,iviktT2) 



and 



Ef nst \ T ,v,k)= ' 



k r Q-i,iw{k t T) 
vQi,i V (k t T) 



K-T Si yiv (k t T) 
k%VTS-lj v (ktT) 



p(r, v, k) 
p(r, v, k) 



(A.21) 



(A.22) 



Once again, in the static limit, j m — 0, and p — 0; thus, — and only the second term in E^ nst ^ survives and 
becomes the previously found E^ stat]l . Notice that the time integration in expressions for potentials looks as retarded, 
T\ > t 2 . This has nothing to do with causal (and the only one meaningful) retardation. This inequality is due to the 
boundary conditions imposed on Ai (to fix the gauge ) when Ai is being rebuilt from Ei, which is simultaneous with 
the sources. The same inequality appears when we shall rebuild the charge dencity p(r) via j m (r)at the previous 
time. 

Now we can move to the fields produced by the currents, and leaving the vanishing effect of p[r — 0) aside. We 
want to present the propagator in its general tensor form which implies that 



[long] 



( Xl )= I d i x 2 Dt n9 \x ll x 2 )r(x 2 ) 



A 



Let us begin with the electric fields produced by the component j n of the current: 



s\,iu{ktt)tdt 



(A.23) 



E)> nst \T U v,k\f) 



(2tt) 3 



n [1 — 8i t iv(ktTi)] / T 2 dT 2 j v (Ti,iy,k) 



n j n (n,v,k) 



~t 2 r T1 

■ — / si, iu (k t t)tdt 
. 1 Jo 



(A.24) 



We see, that the e\^ cancel out the second term in E^ nst ^ , originating, in its turn, from the term with p in Eq. ( A.22] ). 
In this way, we obtain the full form of the 77-component of the longitudinal field, 



E^(T,U,k\f) 



T + TSl til/ (k t T) 



T r 2 dr 2 

3J '[T2,V,k) , 







A^( n ,u,k\f) = 



r 2 dr 2 

w? 



si,iu{ht)tdt 



f'(T 2 ,i/,k) , 



(A.25) 
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where the first term is independent of v and fc, and yields the contact part of the propagator, which (in the coordinate 
representation) reads as 



D^^in^m - W fi - ri) = 



-S{m - - ri) 



(A.26) 



The first line of Eq. (A. 24) clearly illustrates its origin: We started in Eq. (A. 22) with the product 
vs-i y i V {ktT)p{T,v,k). Then, since p(r) is developed dynamically, we expressed p via d^p — > vf \ gaining an ex- 



tra power of v. This allows us to use the relation between two Lommel functions, Eq. (A. 5), and replace (in fact, 
after integrating by parts) v 2 s-\.iu — » 1 — which is equivalent to a straightforward account for the Gauss law. 

The v- and fct-independent unit gives Eq. ( |A.2(i| ).p| 

The second integral term in Eq. ( A. 25 ) can also be Fourier- transformed into the coordinate representation. We 
want to do that here, in order to verify that the contact term is not singled out artificially, and that it is not canceled 
by something hidden in the second term. To compute the integrals from the function S\^ v it can be conveniently 
decomposed in the following way, 

si,i V (x) = Sx t i„(x) - h iv {x) , 

hiu(x) = e '7 /2 ■ l V ' 2 l0 M l ) (x)+H^l{x)] . (A.27) 
2 smn(7w/2) 

The function h^(x) obeys the homogeneous Bessel equation, and thus can describe the field only outside the domain 
of the source influence. In the course of calculations, we use the following integral representation for the Hankcl 
functions, 



- nv/2 H® (fc t r) = — r e T*k t Tcoshe e ^e d9 
n J-oo 



The Lommel function Sn v has a similar representation, 



-x sinh u 



du 



Si 1 iv(x) = x / coshucos^ 
Jo 

which allows one to compute the integral dv exactly, 

/DO 
S^ lv (k t T)e lvr, dv = TrfctTCOshTye^™ 11111 " 1 
-oc 



and from Eq.(A.28) it follows, 

dve l ^h lv {k t r) 



d 



sin[fc t TCOsh6')] 
cosh 2 (6> + 77) 



d 6 



sin[fc t r cosh(0 — 77)] 



cosh 2 9 



(A.28) 



(A.29) 



(A.30) 



(A.31) 



12 One may wonder, why the same type contact term does not show up in other dynamics (and gauges, like A — 0). The 
propagators of these gauges are constructed in such a way that the translation invariance and the possibility of a simple 
momentum representation are preserved. The price for this apparent simplicity is the spurious poles in the propagator without 
a physically motivated prescription to handle these poles. These poles reflect an intrinsic uncertainty in the way one can 
approach the limit of the static field. In order to fix the gauge A = completely, one has to impose some boundary condition 
on the gauge fields at some time t, thus corrupting the translation invariance and gaining additional terms in the propagator, 
which,in fact, are of the same origin as the contact term in the gauge A T = 0. At large t\ and T2, and locally in the coordinate 
rapidity r\ (when the curvature of the hypersurface of the constant r becomes negligible) , the gauge A T — can be approximated 
locally by the gauge A = ||, provided the boundary conditions at r = are released. In this domain, the contribution of 
the contact term is suppressed by the two small curvature factors, g r/rl {Ti)g r,r, (T2) = t^ 2 ^ 2 - Therefore, if a usual scattering 
process between the asymptotic states takes place at large r, wedge dynamics will treat it according to the standard scattering 
theory. 
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Next we may write the full Fourier transforms. From Eq. ( A.30| ), we have 



d 2 k 
(2tt) 3 



. //. /•' 



V 2 



r cosh 7/ 

47T 



t cosh n _ 



, J Q (kr)e- ktTsinh ^dk 
o 



1 



(r 2 + t 2 sinh r/) 1 / 2 



(A.32) 



Starting from Eq. ( A. 31 ), we continue by introducing k z — k t sinh (9 and fco = k t cosh (9 = |k| and changing d 2 kd9 for 
the three-dimensional integration d 3 k. With t = t cosh 77, r = (x, y, r sinh 77), this leads to 



d 2 k 
(2^)3 



dve ivr >h iv {k t T) 



-vi 



d 3 k 



(2tt) 3 |k| 3 



sin fcgt 



d 2 k 
2i(27r) 3 

v 2 r 

V± 0(r 



30 dfc 

T3 

00 



2 1 gikot—ik z z 



4?T 



2 ^ 



r cosh 77 



(r 2 -r 2 ) 



(A.33) 



(r 2 + t 2 sinh 2 rj) 1 / 2 

Adding (|A.32| ) and flA.33p we indeed find that the 7777 component of the longitudinal propagator vanishes at the 
distances r t exceeding r, i.e. outside the light cone of the position of the current which creates the field. Finally, 



^](r,r ?1) fi|i")= / d m df : 



(2tt) 3 



47T 



0(t - r t ) 



f'(T2,V2,r 2 )<, - TS(rj)S(r) 



t 2 sinh rj 



(r 2 + t 2 sinh rj) 1 / 2 



(A.34) 



where, 77 = 771 — 772 and r = r\ — T2- 

The first (contact) term in this formula is indeed very special. It is not limited by the light-cone boundary. The 
second term, does have this boundaries, which are just imposed on the Coulomb-type fields rewritten in terms of 
the natural coordinates of wedge dynamics. It also includes the radiation fields propagating along the light cone 
t = r t . Therefore, only this term can interfere with the radiation fields of the transverse modes. This is clear evidence 
that the cancellation between the contact term and the non-local parts of the propagator is impossible. As it was 
demonstrated in Ref. |2|, the transverse electric field is governed by a usual relativistic wave equation. Integrating 



Eq. (A.34) over r from to n, we recover the potential, and the 7777 component of the propagator, 



4T 9l (Ti,r 2 ;77,rO 



-S(rj)S(f) 



vi 

47T 



tdtQ{t - r t ) 



t cosh 77 

W) 



The remaining components of the propagator A|^™ s ' are 



t cosh 77 

W) 



(A.35) 



(A.36) 



^°" 9l (n,r 2 ;7,,r) = ^|^-r t ) 



Tl r 2 sinh 77 dt f Tl t 2 sinh 77 dt 



R3(t) 



RHt) 



+ tanliT? J 2 S(t - r t )dt\ = (r 2 , n; -77, r) 



(A.37) 



where R{t) — [r 2 +t 2 sinh 2 77] 1 / 2 . The propagator identically vanishes at r t > r, and the derivatives of the step-function 
are confined to the light cone corresponding to the transient radiation which accompaties the creation of the color 
charges. 
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VI. APPENDIX B. SUBLEADING TERMS IN FORWARD SCATTERING OF THE BOOST STATES. 



In the limit of the forward scattering, the general formula (3.15), can be rewritten in the following form, 



Mi, a _n',a' = Sr s ( v i +V2-v[- v' 2 )8(h + k 2 -k[- k' 2 ) 

/ ridri / r 2 * 2 i?W(m;r 1 )^ 2 i ) (™ 1 r 1 )i/«(m' 2 r 2 )^ 2 2 ) (m 2 T 2 ) 
Jo Jo 12- 

+ n [oo]/ . /" ^\ s r n[oo]fn- n- -r rt\ V2 + v Jl 

H 2 — D r,s (n,r 2 ;c,Q) g -g ^„ J (n,T 2 ;C,gj — 2 — 



(B.l) 



where we took the initial transverse momenta Pi = p 2 = 0, and correspondingly, the final state momenta, j} 2 — —p\ = 
q. By its design, the full T-ordercd propagator is a sum of the longitudinal part and two terms originating from 
the transverse electric, V^ TE \ and transverse magnetic, V^™\ modes of the radiation field, 

£>[00] = D [00,long] + jj[QQ](TE) + £) [00](TM)_ 



The r\r\ component of the longitudinal part of the propagator can be read out from the Eq. (A.25) 



4'°"%i,t 2 ;C,£) = — 



si,ic(ltt)tdt 



(B.2) 



where the first term on the right has already been used in Eq. ( 3.20 ) to obtain the main estimate ( 3.27 ). The second 
term yields 



f T d n f T dT 2 f n y iC n 

h= / —J sign(ri-r 2 ) / Si^{q t t)tdt 



T 2 \T 2 



The behavior of the Lommel function in the limit £ — > can be found from the integral representation (A.4) 



Si,io(qtr) = - [1 - cos(/c t rsin0)] d<f> = l- J (qt T ) 



Expanding the Bessel function Jo(qtt) at small qt, and integrating, we arrive at 



2h = 



1 



sin[Cln(T/r )] 



- 1 



C 2 + 16 [ V~ ' T 4 / c 

This term vanishes in the limit of the forward scattering, q% — ► 0. 
The rs component of the longitudinal field propagator, 



D^ n9] (n,r 2 ;C,q) 



J_ qr<h 

2ir qf 



brings in the term 



h = Qt I T i dl ~i I T 2 dT 2 



7^V + cos[Cln(T/r )] 
T 4 



dt 

si,idqtt)— , 



dt 



(B.3) 



(B.4) 



(B.5) 



sign(ri-r 2 ) / Si^(q t t) 



At C — ► and at small qt, it can be represented as the integral 

2I 2 



g 4 T 6 ,1 



y 5 dy / [x l<: + X -^ }{l-x 2 )dx , 

ro/T Jr /Ty 



(B.| 



(B.7) 



(B. 
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which also vanishes in the limit of the forward scattering, q t — > 0. 

The contribution of the components D^" 9 ' and D^° ng ^ into the matrix element (B.l), as well as of all components 
dF W ' of the transverse propagator is estimated exactly in the same way. All these components are defined as the 
integrals, from to rq t , of the functions which are regular at the origin. Therefore, at small qt, all these terms have 
at least one factor q\ and vanish in the limit of the forward scattering. 

The only exception from this scheme is the piece connected with the transverse part D^ TE ^ of the propagator. This 
part is the bilinear form D^ E ^ which has only rs-components and includes the projector S rs — q r q s jq\. Since this 
projector is orthogonal to the vector q t , the contribution of this mode to the matrix element (B.l) identically vanishes. 
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